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ABSTRACT 


Calculating  the  lift,  drag,  and  pitching  moment  coeffi- 
cients for  an  airfoil  from  the  static  pressure  distribution 
obtained  from  wind  tunnel  tests  is  routine  task  when  steady 
flow  is  considered,  but  it  is  much  more  complicated  when  the 
airfoil  is  operating  in  an  unsteady  flow  field,  similar  to 
that  experienced  by  a helicopter  rotor  blade,  produced  by 
an  oscillating  wind  tunnel.  A data  reduction  routine  capable 
of  condensing  the  large  numbers  of  data  associated  with  the 
unsteady  investigation,  as  well  as  a numerical  integration 
algorithm  for  the  unsteady  aerodynamic  coefficients,  were 
developed;  however,  no  unsteady  data  were  collected  due  to 
hardware  failures.  The  ability  of  the  program  was  demonstrated 
on  previously  obtained  steady  and  quasi-steady  data  and  sample 
results  were  presented. 
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I.  INTRODUCTION 


The  primary  objectives  of  this  investigation  were: 

1)  to  report  on  the  data  reduction  techniques  used  to  ob- 
tain values  for  both  steady  and  non-steady  aerodynamic 
coefficients,  resulting  from' the  experimentally  acquired 
static  pressure  distribution  over  the  model  when  the  tunnel 
was  operating  at  a specified  dynamic  pressure,  with  Coanda 
blowing  oscillated  at  various  frequencies  and  amplitudes  of 
momentum  blowing  coefficient;  and  2)  to  present  a sample 
of  results.  Removal  of  the  model  from  the  tunnel  for 
modifications  in  November  1976  precluded  any  actual  process 
ing  of  non-steady  data  consistent  with  the  stated  goals. ^ 
However,  the  data  reduction  techniques  are  applicable 
to  the  pressure  distribution  obtained  with  steady  Coanda 
blowing  and  will  be  demonstrated  on  such  data  obtained  in 
April  1976.  Prior  to  removal  of  the  model,  the  tunnel  was 
operated  without  blowing  in  an  effort  to  generate  a set  of 
non-steady  data  to  be  used  in  tailoring  the  data  reduction 
procedures  to  be  used  downstream  in  the  on-going  project. 

In  that  test,  the  model  experienced  self-excited  harmonic 
pressures  resulting  from  vortices  being  shed  alternately 


^The  modifications,  to  correct  for  surface  delamina- 
tions near  the  upper  jet  nozzle,  include  a steel  upper 
surface  trailing  edge  and  slot  height  adjustment  to  0.015 
inches.  Additional  pressure  taps  were  installed  on  the 
upper  surface  forward  of  the  Coanda  slot,  and  leaking 
tubes  were  repaired  while  the  modifications  were  being 
implemented. 
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from  the  upper  and  lower  surfaces.  The  manually  recorded 
data  consisted  of  the  RMS  and  mean  values  of  the  unsteady 


pressure,  and  the  phase  angle  in  respect  to  a reference  for 
each  pressure  orifice.  A more  complete  description  of  the 
data  acquisition  system  is  given  in  Section  II  B. 


II.  EQUIPMENT  AND  INSTRUMENTATION 


A.  AIRFOIL 

The  airfoil  used  in  this  study  was  obtained  from  the 
Lockheed-California  Company  after  completion  of  their 
design  feasibility  study  on  circulation  control  rotors 
(Ref.l).  It  had  an  elliptical  cross-section  with  a chord 
of  10.215  inches,  thickness  ratio  of  0.20,  and  camber  of 
3.3%  (Figure  1).  A slot  to  facilitate  Coanda  blowing  was 
machined  in  the  upper  surface  at  chord  position  X/C=0.95, 
with  the  slot  height  set  at  0.010  inches.  Coanda  air  was 
supplied  to  the  plenum  cavity  at  both  ends  of  the  section 
to  insure  uniform  flow  along  the  span.  Fifty-three 
pressure  taps  at  mid-span  with  a chord-wise  orientation  were 
monitored  through  two  scanivalve  transducers  for  an  elec- 
trical representation  of  the  pressure  distribution  over  the 
airfoil.  The  pressure  tap  locations  are  listed  in  Table  I. 

The  2.0  X 2.0  foot  unsteady  flow  wind  tunnel  in  the 
Department  of  Aeronautics  Laboratories  at  the  Naval  Post- 
graduate School  was  used  as  it  is  capable  of  simulating  the 
velocity  flov-'  field  experienced  by  blade  elements  of  a 
helicopter  in  forward  flight.  The  model  was  mounted  at  an 
angle  of  attack  of  -5  degrees,  which  closely  corresponds  to 
the  zero-lift  angle  of  attack  for  the  data  reported  herein. 


FIEURC  1 

RIRFDIL  CRD55-5EaiDN  WITH  THP  LDCRTIQN5 


TABLE  I 


HlRFDiL  PRESSURE  TAP  LOCATIONS 


TAP  * 

X 

Y 

TAP  * 

X 

Y 

1 

0.000 

0 . 000 

20 

10.20E 

“0.020 

2 

0.003 

0.0S0 

29 

10. 1E2 

“0.040 

3 

0.0SS 

0.  13S 

30 

10. 1 10 

“0. 152 

H 

0.  123 

0.220 

31 

10.051 

“0.255 

E 

0.21  1 

0.30S 

32 

9.940 

“0.372 

E 

0.327 

0.H00 

33 

9.000 

“0.400 

7 

0.S2S 

0.S1B 

34 

9.E13 

“0.443 

B 

o.gsM 

0.710 

35 

9.1  14 

“0.525 

3 

1 .HHl 

0.B75 

3E 

0.507 

“0.590 

10 

1 .931 

1 .004 

37 

7.950 

-0.E3B 

1 ( 

2.H3S 

1 .0EE 

30 

7.570 

“0.E5B 

12 

2.BS1 

1 . 102 

39 

7.059 

“0 . 570 

13 

3.95B 

1 .315 

40 

E.033 

“0.700 

IH 

S.0BE 

1 .3E0 

41 

5.  1 IE 

“0.700 

IS 

B.09B 

1 .343 

42 

4.019 

“0.E35 

IE 

7.132 

1 .230 

43 

2.902 

“0 . EES 

17 

7.E27 

1 . 1E2 

44 

2.49E 

“0.E45 

IB 

B.019 

1 .000 

45 

1 .992 

“0.E10 

19 

B.E3S 

0.950 

4E 

1 .470 

“0 . 5E0 

20 

47 

0.9B4 

“0.490 

21 

40 

0.521 

“0 . 300 

22 

43 

0.349 

“0.320 

23 

10.070 

0.340 

50 

0.230 

“0.2E0 

2H 

10. IHO 

0.220 

51 

0.  1 15 

“0. 170 

2S 

10. 191 

0.  147 

52 

0.034 

“0.090 

2E 

27 

10. 199 
10.21S 

0.073 
0 . 000 

S3 

0.010 

•- 0.052 

X IN 

INCHES  AFT 

OF  LEADINE 

EDEE 

Y IN  INCHES  FROM  CHORD  LINE 
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B.  INSTRUMENTATION 


The  data  acquisition  system  evolved  from  a classic 
manometer  board  arrangement  with  readings  hand  recorded 
(Ref.  2)  and  reduced  to  the  present  system,  which  employed 
electrical  read-outs  with  an  existing  Digitec  system  with 
customized  signal  conditioning.  The  digital  multiplexing 
circuitry  in  the  Digitec  system  was  interconnected  through 
the  scanivalve  solenoid  controllers  for  automatic  scani- 
valve  advancing  as  the  output  of  the  digital  voltmeter  was 
sequentially  printed  on  a paper  tape.  Pressure  sensing  was 
accomplished  by  connecting  each  of  the  pressure  taps  to 
one  of  the  two  scanivalves  (Figure  2)  through  tubing  of 
uniform  diameter  and  length.  This  was  done  to  insure  a 
constant  transfer  function  for  all  pressure  taps,  when 
measuring  unsteady  pressures.  Calibration  tests  using  the 
method  reported  by  Johnson  (Ref. 3)  showed  an  acceptable 
variation  in  dynamic  gain  in  the  frequency  range  from  0 to 
300  Hz  (Figure  3) , which  included  the  range  of  practical 
interest. 

For  the  unsteady  aerodynamics  analysis  (C^=0) , the 
voltage  representation  of  the  surface  static  pressures  was 
of  the  form 

6x=ax+bx(t).  (2.1) 

Full  analog  measuring  techniques  were  employed  in  measuring 
the  mean  voltage,  RMS  voltage,  and  phase  relationship  to  a 
reference  for  each  pressure  tap.  The  mean  voltage  was  ob- 
tained by  processing  the  signal  through  a low-pass  filter 
circuit  (Figure  4)  with  a two  second  time  constant  which 
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corresponds  to  a corner  frequency  of  0.08  Hz  on  a Bode  plot. 
Therefore,  the  influence  of  the  unsteady  voltage  was  essen- 
tially negated  and  the  output  readings  represented  the 
steady  components  of  the  signal. 

C.  ADVANCED  DATA  ACQUISITION  SYSTEM 

To  meet  the  primary  objective  mentioned  in  Section  I, 
a more  sophisticated  data  acquisition  and  reduction  system 

2 

will  be  required  and  is  now  being  constructed  by  Englehardt 
at  the  Naval  Postgraduate  School  in  the  Department  of  Aero- 
nautics. The  system,  based  on  the  Intel  MDS-80  micro- 
processor unit,  will  control  the  logic  of  the  experiment, 
collect  and  store  data,  and  handle  the  entire  data  reduction 
process.  The  algorithm  for  data  reduction  discussed  in 
Section  V was  proven  sufficiently  simple  and  capable  of 
performing  all  required  data  reduction  and  numerical  integra- 
tion. 


“Master  of  Science  thesis  bv  Englehardt  with  a complete 
description  of  the  entire  microprocessor  system  will  be 
published  in  June  1977. 


C^=20  >ifd 


Figure  4 

Low-Pass  Filter  Circuit 
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III.  DATA  REDUCTION 


A.  STEADY  DATA  PROCESSING 

The  aerodynamic  coefficients  are  defined  by  integral 
equations  but  must  be  approximated  by  numerical  integra- 
tions, since  data  are  available  only  at  a finite  number 
of  pressure  tap  locations.  The  approximate  steady  normal 
force,  chord  force,  and  pitching  moment  coefficients  are 


1.0 


Y/C  (max) 


-f 


Y/C  (min) 
Y/C  (max) 


C3.1) 


'm(TE)  - J + 

-Y/C  (min)  ^ ^ 


1.0 


(C  -C  )(l-^/.)d(^/ 
J - Pi  Pu  ^ 


C) 


Then,  with  an  angle  of  attack  correction  for  the  force 

coefficients  and  a moment  transfer  to  the  quarter- chord 

position,  the  aerodynamic  coefficients  are 

C,  =C„cosoi-C„sina 
In  c 

C , = C„sina+C^COSa  (3.2) 

one 

‘^m(c/4)^^m(TE)’°’'^^^n- 

The  pressure  data  from  the  April  1976  tests  and  the 
November  1976  mean  pressure  data  were  transferred  to  a 
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Hewlett-Packard  9830  tape  cassette  and  processed  using  the 
program  listed  in  Appendix  A,  which  employs  the  equivalent 
numerical  relationship  to  the  equations  stated  above. 


B.  UNSTEADY  DATA  PROCESSING 


Before  the  unsteady  data  were  suitable  for  numerical 
integration,  a dynamic  gain  correction  was  necessary  to 
account  for  the  transfer  function  discussed  in  Section  II. 

The  frequency  of  the  self-induced  oscillations  was  constant 
at  262  Hz  when  the  tunnel  airspeed  was  approximately  100 
ft/s.  From  Figure  3,  the  corresponding  gain  correction, 

|G^|,  was  approximately  0.5,  which  implied  that  the  amplitude 
of  the  surface  static  pressure  was  attenuated  by  a factor 
of  one-half  at  the  transducer.  All  phase  angle  values  were 
relative  to  tube  number  7,  which  was  arbitrarily  chosen  as 
a clock  reference.  Phase  angle  errors  resulting  from  the 
two  different  scanivalve  circuits  were  faired  out  in  the 
analysis  (Figure  5). 

With  the  above  corrections,  the  numerical  integration 
of  the  unsteady  static  pressure  distribution  was  carried 
out  in  a manner  similar  to  the  steady  pressure  integrations. 
For  the  unsteady  normal  force  coefficient 


(C  sin(!>-C  sin<t>)d(^/p) 

Pi  Pu  ^ 


COS(Ji 


1.0 

0 


cos(ji-C 

P 


coS(|i)d(^/„) 
u ^ 


(3.3) 
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Then 


And  = tan"'  ( <:„^>sin») 

( (C^  > cos4i) 

/ 2\^ 

where  \C^  / is  the  RMS  value  of  the  normal  force  coefficient 
and  (j>^  is  the  phase  angle  relative  to  tube  number  7.  The 
chord  force  and  pitching  moment  coefficients  were  calculated 
similarly. 


2 H 

(C^  ) sin(()  = 


Y/C (max) 


- Y/C (min) 


(C  sin<|i-C  sin<())d(  / ) 
Pf  Pr 


Y/C (max) 


/ 2v^  f Y 

<C  ) cos<j)  = (C  coS(j)-C  cos()))d(  /„) 

^ J Pf  P-r  ^ 


•Y/C (min) 


(<C^  > sin<J)) 
(<C^  > cos<(i) 


^ Y/C (max) 

(Cj^^)  sin<i>  = J (C  sin*-C  sin*)  (''^/^) d('^/^)  -P 

-Y/C (min)  ^ ^ 


i. . u 

I (Cp  sin*-Cp  sin*)(l-Vc)d(Vc) 

^ n 1 U 


Y/C (max) 


> cos^  = 


te 


Y/C(min) 


(Cp^cos<^-Cp  cos$)  C^/^)d(^/^)  + 


1 

(C  cos^-C  cos4.)(l-^/  )d(^/„) 
Pi  Pu  L L 


h 


(<Cm  >sin^) 
«Cj^  > cosij>) 


The  unsteady  aerodynamic  coefficients  are  calculated  as  in 
Eqs.  3.2,  with  the  aid  of  elementary  trigonometric  identities 
to  account  for  the  phase  angles.  The  program  listed  in 
Appendix  A has  the  option  for  unsteady  numerical  integration 
incorporated. 

It  is  intended  that  the  C^=0  oscillating  pressure 
measurements  be  repeated  when  the  improved  data  acquisition 
system  becomes  available.  The  tests  will  be  conducted  at 
several  tunnel  dynamic  pressures  to  determine  the  Reynolds 
number  dependence  of  the  data.  The  previous  measurements 
were  at  approximately  R^=0.54  x 10^. 


IV.  RESULTS 
A.  STEADY  DATA  RESULTS 

The  validity  of  the  concept  of  circulation  control  was 
demonstrated  by  Englar  (Refs.  4 and  5)  and  is  confirmed 
here  with  the  steady  momentum  blowing  coefficient  results. 
The  resulting  aerodynamic  coefficients  are  summarized  in 
Table  II.  The  ability  to  produce  relatively  large  in- 
creases in  lift  coefficient  with  only  a small  increase  in 
blowing  coefficient,  shown  graphically  in  Figures  6a,  b, 
and  c,  is  one  of  the  desirable  characteristics  of  the 
system.  Probably  the  best  measure  of  the  ef fectiveness  is 
the  aerodynamic  lift  amplification  ratio,  3C^/8C^,  which 
was  computed  to  be  approximately  15.3,  the  low  value  being 
attributed  to  the  geometric  deficiencies  of  the  trailing 
edge  slot  as  alluded  to  in  Section  I.  However,  the  results 
are  encouraging. 

Close  scrutiny  of  the  coefficient  of  pressure  informa- 
tion plotted  against  X/C  (Figures  7-12}  reveals  a pressure 
distribution  over  the  model  such  that  the  center  of  pressure 
was  at  the  mid-chord  position,  which  is  in  agreement  with 
the  theory.  The  center  of  pressure  position  is  also 

verified  by  the  graph  of  C , Vs.  C . 

° ^ m(c/4j  u 

rX-,  C (3C  /3C  } (-3.5} 

(tt}  = t - m'  u''  = . 25  - •'  = . 48 

^ C.P.  ^ (3Cj/3C^}  (TSTT) 
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TABLE  II 


AERQDTOAMIC  COEFFICIENT  3UI>:MARY 
oTEADY  AEROBYNAMIC  C0EFPICIEKT3 
Angle  of  Attack=-5  degrees  U=100  ft/s 


Run  No. 

C/A. 

^1 

°d 

*^m(c/4) 

42301 

0 

-0.0942 

0.0341 

-0.1065 

42302 

0.0855 

1.0332 

0.1083 

-0.3907 

42303 

0.0640 

0.8146 

0.0678 

-0.3365 

42304 

0.0448 

0.5774 

0.0440 

-0.2776 

42305 

0.0272 

0.3459 

0.0327 

-0.2175 

432O6 

0.0133 

0.1100 

0.0367 

-0.1664 

UNSTEADY  AERODYTJAMIC  COEFFICIENTS 

C^ =-0.0765  C^(RM3)=2.1^04  i>  =225 

C^=  0.0451  C^(HiM  3)  =0.1629  = 59 

G =-0.1063  C„(  RMS)  =0.5817  ^=65 

in  in 


degrees 

degrees 

degrees 
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FIGURE  BC 
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CP  V5.  X/C 


3 


Also,  the  low  drag  qualities  of  the  airfoil  are  favorable 
when  full  thrust  recovery  of  the  momentum  blowing  coefficient 
is  assumed.  Pressure  coefficient  data  are  listed  in 
Appendix  B. 


B.  UNSTEADY  DATA  RESULTS 


The  unsteady  aerodynamics  of  the  elliptically  shaped 
airfoil  have  been  given  a rather  cavalier  treatment  in 
the  past  but  their  contribution  to  the  total  forces  and 
moments  must  be  realized.  The  implications  of  the  unsteady 
aerodynamics  is  not  generally  well  understood;  therefore, 
a simple  example  is  presented  here  to  emphasize  their 

I 

importance  before  discussion  of  the  unsteady  results. 

Consider  an  airfoil  at  zero  angle  of  attack  with  the 
t RMS  pressure  distribution  equal  on  the  upper  and  lower 

j 

surfaces  (Figure  13a)  but  with  unequal  relative  phase 
angles  (Figure  13b)  taken  from  the  quarter-chord  and  the 


three-quarter-chord  positions  from  the  actual  unsteady  data 

/ 2\'^ 

(Figure  5).  The  unsteady  normal  force  coefficient,  / , 
is  calculated  from  the  numerical  approximation  to  Eq.  3.3, 
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Similarly,  for  the  moment  coefficient  (referenced  to  the 

trailing  edge) : 

- % 

sin<^  = 0.5((sin(210°)-sin(35°))0.75+(sin(230°)-sin(55°))0.25) 


2. 


siniji  = -0.6 


cos  * 0.5((cos(210°)-cos(35°))0.75+(cos(230°)-cos(55°))0.25) 


> cos  = -0.78 


= "^0.6^  + 0.78^ 


-0.6 

-OTTF 


<C  =0.98 
' m ' 


4.  = 217' 


Clearly,  the  force  and  moment  coefficients  are  of  significant 
magnitude  even  if  a uniform  distribution  (RMS)  of  pressure 
coefficient  were  assumed,  because  of  the  importance  of  con- 
sidering relative  phase  relationships.  The  dynamic  behavior 
of  the  unsteady  aerodynamics  of  the  circulation  control 
airfoil,  particularly  the  lift  amplification  ratio,  in  the 
unsteady  environment  caused  by  cyclic  blowing  will  be  in- 
vestigated after  the  improved  data  acquisition  system  is 
operational . 

As  discussed  in  Section  I,  measurements  were  taken  in 
November  1976,  with  the  airfoil  experiencing  harmonic 
pressure  variations.  The  unsteady  data  are  listed  in 
Appendix  C.  The  experimental  data  were  numerically  inte- 
grated, similar  to  the  previous  example,  for  the  lift,  drag, 
and  pitching  moment  coefficients  (Figures  14  and  15) . It 
should  be  recognised  that  the  phase  shift  between  the  upper 
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1 


and  lower  surfaces  is  approximately  180  degrees,  the  result 
being  the  relatively  large  amplitude  of  the  aerodynamic 
coefficients.  Highlights  of  the  unsteady  results,  summarized 
in  Table  II,  were  that: 

1.  The  time  varying  pressures  over  the  model  were 
predominantly  harmonic,  with  a non-dimensional  frequency 
(Strouhal  no.)  of 


2.  The  RMS  value  of  the  normal  force  coefficient  was 
approximately  2.14,  which  was  an  order  of  magnitude  greater 
than  the  RMS  drag  coefficient. 

3.  The  test  Reynolds  no.  was  approximately 

R^  = 0.54  X 10^. 
e 

4.  The  location  on  the  airfoil  for  the  center  of 
pressure  due  to  unsteady  loadings  was  approximately  0.491 
chord  and  the  RMS  value  of  pitching  moment  at  this  reference 
point  was 

2 ^ 

<C  = 0.1694 
' m ' 

with  a phase  angle  of  90  degrees  lag  relative  to  the  normal 
force.  This  term  may  be  considered  as  an  apparent  mass 
type  effect. 
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V.  DATA  REDUCTION  FOR  THE  ADVANCED  DATA  REDUCTION  SYSTEM 


The  advanced  data  acquisition  system  described  in 
Section  II  C will  allow  a more  comprehensive  investigation 
o£  the  unsteady  aerodynamics  of  the  model,  but  the  data 
reduction  techniques  must  be  tailored  to  handle  the  large 
numbers  of  data  to  be  collected.  Oscillation  of  the  plenum 
pressure  as  described  by  Bauman  (Reference  6)  induces 
harmonic  blowing  from  the  trailing  edge  slot  at  a specified 
frequency,  It  is  assumed  that  the  unsteady  static 

pressure  at  each  station  will  have  a steady  value  plus  a 
periodic  value  with  frequency  superimposed  due  to  the 
oscillating  momentum  blowing  coefficient. 

Since  the  data  are  assumed  periodic  in  form,  the  un- 
steady pressure  at  any  point  on  the  airfoil  may  be  defined 
by  a Fourier  series  (Ref.  7)  as  follows: 

P(x,t)=A^(x)+^|^  cos(nW^t)+Bj^(x)  sin(nW^t)]  (5.1) 

where  is  the  frequency  of  the  oscillation.  The  standard 
formulae  for  obtaining  the  Fourier  coefficients  for  the 
infinite  series  are: 

T 


A,(x)  = f P(x)dx 

0 


2T 

T 


Aj^(x)  = ^ 


P(x)cos (nx)dx 


(5.2) 


-T 


1 f T 

Bn(x)  = Y P (x) sin(nx) dx 
■'  -T 
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Because  the  experimental  data  are  in  the  form  of  a 
discritization  of  a truncated  periodic  signal,  the  above 
equations  are  slightly  altered  for  estimating  the  Fourier 
coefficients.  The  modified  form  is: 


X^(x) 


V 1 ^ 

= 4 z:  p.  (X) 

o i-1  ^ 

2 ^ 

= w r P.  (x) cos  (nW,  t)  ; n=l 
^ i=l  ^ 1 


=1-2:  P.  (x)sin(nw  t) 

n ^ i=i  ^ 1 


,n  (5.3) 


where  k is  the  number  of  samples  to  be  analyzed  and  n 
is  the  harmonic  number.  It  is  recognized  that  the  physical 
meaning  of  A^(x)  is  the  steady  or  mean  value  of  the  oscillat- 
ing pressure.  That  is,  A^(x)  = ^p^(x))  where 


Lim 

T-V  00 


1 

It 


J 

P (x) dx 
-T 


(5.4) 


By  subtracting  out  the  mean  pressure  value,  the  signal  becomes 
biased  to  a pure  oscillating  pressure  with  a mean  value  equal 
to  zero.  The  terms  A^^(x)  and  Bj^(x)  are  equivalent  to  the 
cosine  and  sine  components,  respectively,  discussed  in 
Section  III  B.  Figure  16  illustrates  the  reconstruction  of 
a discritized  signal  utilizing  this  algorithm  programmed  on 
the  Hewlett-Packard  9830  calculator  (program  listing  in 
Appendix  A).  The  error  was  minimized  by  insuring  that  the 
number  of  samples  analyzed  was  an  integer  multiple  of  the 
frequency;  however,  the  truncation  error  was  never  greater 
than  2%. 
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PLOT  or  5INCX)+.B«51NC2*X) 


3anss3yd  x<iy3iSNn 
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Having  developed  an  effective  method  for  handling  the 
large  numbers  of  data,  calculation  of  the  aerodynamic  force 
and  moment  coefficients  can  be  accomplished.  Before  the 
actual  processing  of  data  begins,  the  sample  pressures  are 
converted  to  non-dimensional  coefficient  form  to  expedite 
calculations.  The  unsteady  pressure  coefficients  are 


i = 1,2,3,... ,53  (5.5) 

where  P and  H are  steady  pressure  values  from  a pitotstatic 

tube  mounted  on  the  floor  of  the  wind  tunnel.  Then  the 

oscillating  signal  at  any  pressure  tap  is; 

C (x,t)  = l A cos(nw,t)+B  sin(nw.t) 
p ’ n I n 1 


Cp(x,t)  = 2 sin(nw^t  + 4>j^) 


(5.6) 


where 


T = "mA  ^+B  ^ ^ 

n in  n ’ ^n 


n 


The  unsteady  C^(x,t)  are  related  by  the  phase  angles  (}i  , 
however,  the  phase  angles  are  meaningless  unless  they  are 
measured  relative  to  a "clock"  or  reference  signal.  There- 
fore, it  is  necessary  to  designate  one  pressure  tap  to  be 
used  as  the  reference.  Then,  the  relative  phase  angle  is 

(5.7) 

^rel  ^ref  ^ ^ 

For  simplicity,  shall  hereafter  be  interpreted  as  the 
relative  phase  angle. 

With  the  steady  pressure  coefficient  removed  from  the 


oscillating  signal,  separate  aerodynamic  coefficients  can 
be  calculated  for  the  steady  and  unsteady  signals.  The 


I 


f 


steady  force  and  moment  coefficients  are  found  exactly  as 
in  Section  III  A.  The  unsteady  aerodynamic  coefficients  can 
be  defined  in  a Fourier  series  form  as: 

C^(t)  = Z sin(kw^t  + (()j^) 

C^(t)  = Z sin(kwj^t+(|)j^)  (5.8) 

^m^^^  = ^ ^m^  sin(kw^t+<|.j^) 

The  magnitudes  and  phase  angles  of  the  coefficients  can  be 

determined  thusly: 

for  the  first  harmonic 

sin(w.t  + (}i)  = C sin(|)Cos  (w,  t) +C  cos4>sin(w.  t) 
n^^  1 n.  I n,  1 

where 

C sin({)  = 

^1 


1 

I |Cp^(f)sin»j(|)-Cp^(|)siM„(|)]  d(|) 


C cos((i 
^1 


(|)cos^^(|)  - 


'n  (7)C0S4> 

Pu  ^ 


u'|)] 


) (5.9) 


then 

- ■^(Cn^cos.i)^  . - Tan 

Similar  calculations  are  done  to  obtain  values  for  the  higher 
harmonics  (C^^  ,C^  , * * " ) as  well  as  the  chord  force  and 

pitching  moment  coefficients  using  the  higher  harmonics  of 
pressure  coefficient  information  as  obtained  during  the 
initial  data  reduction  previously  described. 


-1 

c_  cos^. 


n, 
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VI.  CONCLUSIONS 


As  mentioned  in  Section  I,  the  two-fold  objective  of 
this  paper  was  to  develop  a computer  program  for  unsteady 
aerodynamic  analysis  and  to  demonstrate  its  capabilities 
on  experimentally  generated  unsteady  data.  The  objective 
was  partially  satisfied;  a computer  program,  well  within 
the  limitations  of  mini-calculators  such  as  the  Hewlett- 
Packard  9830,  was  developed  which  employed  a straight- 
forward numerical  integration  of  pressure  coefficients  to 
obtain  values  for  the  aerodynamic  coefficients.  The  pro- 
gram can  be  modified  for  unsteady  data  by  properly  taking 
into  account  the  cosine  and  sine  components  of  the  pressure 
signal  relative  to  a fixed  harmonic  (clock)  signal.  Also, 
the  extraction  of  the  Fourier  coefficients  of  the  truncated, 
time-varying  signal  appears  to  be  within  the  capabilities 
of  the  updated  data  acquisition  system. 
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APPENDIX  A 


COMPUTER  PROGRAM  LISTINGS  AND  OUTPUT 


1 REM  PROGRAM  FOR  DECOMPOSITION  OF  A FUNCTION 

2 REM  INTO  FOURIER  COEFFICIENTS 
10  DIM  EC256: 

30  W1=(20*PI) 

35  DISP  "DIVISIONS  PER  CYCLE?"; 

36  INPUT  D 

42  P=2*PI/W1 

43  T=P/D 

50  REM  COMPUTE  SAMPLE  FUNCTION 
60  FOR  1=1  TO  220 
65  X=W1*(I-1)*T 

70  ECI3=SIN(X)+0.5*SIN(2*X)+0.3*SIN(3*X) 

80  NEXT  I 
100  PRINT 

105  PRINT  "THE  DRIVING  FREQUENCY  IS  10  HZ" 

106  PRINT  D;  "DIVISIONS  PER  CYCLE" 

107  PRINT 

110  PRINT  " COEFFICIENTS" 

111  PRINT 

113  PRINT  TAB2"COS"TAB15"SIN" 

114  FOR  J=1  TO  3 

115  PRINT 

117  PRINT  TAB2"A";J,  "B";J,  "PHASE  ANGLE" 

119  A1=B1=0 

120  FOR  1=1  TO  220 

130  Bl=(Bl  + ECn*SIN(J*Wl*(I-l)*T)) 

140  Al= (Al+ECn  *COS(J*Wl* (I- 1) *T) ) 

150  NEXT  I 

159  P=ATN(B1/A1) 

160  WRITE  (15,170HA1*2/220)  , (Bl*2/220)  , (P*180/PI) 
170  FORMAT  F6 . 3 , 7X,F6 . 3 , 13X, F5 . 1 

175  NEXT  J 
180  END 


THE  DRIVING  FREQUENCY  IS  10  HZ 
10  DIVISIONS  PER  CYCLE 


COEFFICIENTS 


COS 

SIN 

A 1 

B 1 

PHASE  ANGLE 

0.000 

1.000 

90.0 

A 2 

B 2 

PHASE  ANGLE 

0.000 

0.500 

90.0 

A 3 

B 3 

PHASE  ANGLE 

0.000 

0.300 

90.0 

1 REM  COMPUTE  CL, CD, CM- -STEADY  AND  UNSTEADY 

2 REM  FROM  CP  DATA  ON  TAPE  CASSETTE 
10  COM  R1,DC80J  ,PC80J  ,CC60] 

20  DIM  AC1063 
30  PRINT 

40  LOAD  DATA  4, A 
60  C9=0 
70  X9=0.95 
80  Y9=0.048 

90  DISP  "CP  DATA  FILE  NO.  IS"; 

100  INPUT  Z 

110  LOAD  DATA  Z 

111  DC223=  0.4 

112  DC36J=0.539 

113  Gl=l 

115  PRINT  "INPUT  OPTION" 

116  PRINT  "1  FOR  STEADY  DATA" 

117  PRINT  "2  FOR  NON-STEADY  DATA" 

118  INPUT  G3 

119  GOTO  G5  OF  130,212 

120  GOTO  130 
130  PRINT 

140  PRINT  "RUN  NO.  IS"R1 

150  PRINT  TAB7"UPPER"TAB24"LOWER" 

160  PRINT  TAB2"N"TAB9"CP"TAB19"N"TAB26"CP" 

170  FORMAT  F3 . 0 , 2X ,F7 . 3 , 5X,F3 . 0 , 2X , F7 . 3 
180  FOR  N=1  TO  27 
190  M=N+26 

200  WRITE  (15,170)N,C[Nl  ,M,C[M] 

201  NEXT  N 

202  FIXED  3 

203  PRINT  "CP1="C9"AT  X/C="X9",  Y/C="Y9 

204  STANDARD 

2C5  DISP  "CORRECT  THE  CP  VALUES" 

206  STOP 
210  GOTO  270 

212  DISP  "INPUT  DYNAMIC  GAIN"; 

213  INPUT  G2 

214  FIXED  3 

216  PRINT  "RUN  NO.  IS"R1 

217  C9=CL2I] +0.33*(C[22J -CL21]  ) 

218  PRINT  "TAP"TAB7"CP(MEAN)"TAB17"CP(RMS)"TAB27"PHI" 
220  FORMAT  F3 . 0 , 3X, F7 . 3 , 3X, F7 . 3 , 3X , F4 . 0 

222  FOR  I»1  TO  53 

224  WRITE  (15, 220)1, C[IJ,D[IJ, PCI] 

226  NEXT  I 

227  PRINT  "C9="C9 

228  DISP  "CORRECT  THE  CP  VALUES" 

230  STOP 

235  GOTO  270 

240  FOR  1=1  TO  53 

242  C[I]=G2*D[I]*SIN(P[IJ/57.296) 

244  NEXT  I 

245  C9=(D[21]  + 0. 33*(D[2  2]  -D[21]  )) 

246  Gl=2 

247  GOTO  270 
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250  FOR  I»1  TO  53 

252  C[IJ-G2*D[I]*COS(P[I]/57.296) 

254  NEXT  I 

256  C9»(D[21J+0.33»(D[22]-DC21])) 

257  Gl=3 

258  GOTO  270 
270  L1»0 
272  L2=0 
274  P1=0 
276  P2»0 

310  FOR  1=1  TO  26 
320  IF  1=21  THEN  490 
330  C1  = -C[IJ 
340  C2  = -C[I  + 1J 
350  X1=A[IJ 
360  X2=A|:i  + ]J 
370  GOSUB  1120 
380  L1=L1+N1 
390  P1=P1+M1 
400  Cl  = C[I  + 27] 

410  C2=CCI  + 263 
420  Xl=ACI  + 27] 

430  X2=ACI  + 261 
440  GOSUB  1120 
450  L2=L2-*-Nl 
460  P2=P2+M1 
470  NEXT  I 
480  GOTO  640 
490  C1  = -C[21] 

500  C2=-C9 
510  X1=A[21] 

520  X2=X9 
530  GOSUB  1120 
540  L1=L1+N1 
550  P1=P1+M1 
560  Cl=-C9 
570  C2=-C[22] 

580  X1=X9 
590  X2=A  22 
600  GOSUB  1120 
610  L1=L1+N1 
620  P1=P1+M1 
630  GOTO  400 
640  D1=0 
650  D2=0 
660  P3=0 
670  P4=0 

680  FOR  1=1  TO  13 
690  Cl+CCU 
700  C2=CCI  + 1] 

710  Y1=ACI  + 531 
720  Y2=ACI  + 541 
730  GOSUB  1160 
740  D1=D1+N1 
750  P3=P3+M1 
760  NEXT  I 
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770  FOR  1=40  TO  52 
780  C1=C[ID 
790  C2=C[I  + i;] 

800  Yl=A[I+53: 

810  Y2=ACI+543 

820  GOSUB  1160 

830  D1=D1+N1 

840  P3=P3+M1 

850  NEXT  I 

860  FOR  1=14  TO  39 

870  IF  1=21  THEN  970 

880  ci=-cci+i: 

890  C2=C[i:i 
900  Yl=A[I  + 54‘] 

910  Y2=ACI+53D 
920  GOSUB  1160 
930  D2=D2+N1 
940  P4=P4+M1 
950  NEXT  I 
960  GOTO  1200 
970  Cl=-C9 
980  C2=C  21 
990  Y1=Y9 
1000  Y2=A  74 
1010  GOSUB  1160 
1020  D2=D2+N1 
1030  P4=P4+M1 
1040  C1=-CC22] 

1050  C2=-C9 
1060  Yl=AC75n 
1070  Y2=Y9 
1080  GOSUB  1160 
1090  D2=D2+N1 
1100  P4=P4+M1 
1110  GOTO  950 

1120  M1=C1*(X2-X1)*(1-0.5*(X2+X1) 

1130  M1=M1+(C2-C1)*(X2-X1)*(0.5- (X2/3)  - (Xl/6)  ) 
1140  N1=0. 5*(C2+C1)*(X2-X1) 

1150  RETURN 

1160  N1=0.5*(C2+C1)*(Y2-Y1) 

1170  M1=0.5*C1*(Y2-Y1)*(Y2+Y1) 

1180  Ml=Ml+(C2-Cl)*(Y2-Yl)*((Yl/6)+(Y2/3)) 

1190  RETURN 
1200  PRINT 

1205  GOTO  G1  OF  1210,1500,1600 
1210  PRINT  "RUN  NO.  IS"R1 
1220  PRINT 
1225  FIXED  4 

1230  PRINT  "COEFFICIENT  BUILDUP  IS" 

1250  PRINT  TAB2"CNU="L1 
1260  PRINT  TAB2"CNL="L2 
1270  PRINT  TAB10"TOTAL  CN="(LI+L2) 

1280  PRINT  TAB2"CCF="D1 
1290  PRINT  TAB2"CCR="D2 
1300  PRINT  TAB10"TOTAL  CC="(D1+D2) 
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1310  PRINT  TAB2"CMU="P1 

1320  PRINT  TAB2"CML="P2  * 

1330  PRINT  TAB2"CMF="P3 
1340  5jaNT  TAB2"CMR="P4 

1350  PRINT  TAB10"TOTAL  CM=" (P1+P2+P3+P4) "AT  T.E. 
1360  M1=(P1+P2+P3+P4) 

1370  N1=(L1+L2) 

1380  C1=(D1+D2) 

1390  M5=M1-(0.75*N1) 

1400  DISP  "INPUT  ANGLE  OF  ATTACK  (DECS.)"; 

1410  INPUT  A2 

1415  A2=0.01745*A2 

1416  A9=A2 

1420  N5=N1*C0SCA2)-C1*SIN(A2) 

1430  C5=C1*C0S(A2)+N1*SIN(A2) 

1440  PRINT 
1445  C2=A2*57.296 

1450  PRINT  TAB  2"CL (MEAN) ="N5 ;TAB21"CDCMEAN) ="C5 
1460  PRINT  TAB2"ALPHA  ="C2"DEGS" 

1470  STANDARD 

1480  GOTO  G3  OF  90,240,250 

1490  END 

1500  N3=L1+L2 

1502  C3=D1+D2 

1504  M3+P1+P2+P3+P4 

1505  FIXED  4 

1506  PRINT  "CN(RMS)*SIN(PHI)="N3 
1508  PRINT  "CC(RMS)*SIN(PHI)="C3 

1510  PRINT  "CM(RMS)*SIN(PHI)="M3 

1511  PRINT 
1595  GOTO  250 
1600  N4=L1+L2 
1602  C4=D1+D2 

1604  M4=P1+P2+P3+P4 

1605  FIXED  4 

1606  PRINT  "CN(RMS)*COS(PHI)="N4 
1608  PRINT  "CC(RMS)*COS(PHI)="C4 

1610  PRINT  "CM(RMS)*COS(PHI)="M4 

1611  PRINT 

1612  Nl=SQR(N3+2+N4+2) 

1613  P1=ATN(N3/N4) 

1614  X=N3 

1615  Y=N4 

1616  F=P1 

1617  GOSUB  1900 

1618  P1=F 

1620  PRINT  "CN(RMS)="N1,"PHI="P1*57.296 
1630  Cl=SQR(C3+2+C4+2) 

1640  P2=ATN(C3/C4) 

1641  X=C3 

1642  Y=C4 

1643  F*P2 

1644  GOSUB  1900 

1645  P2-F 

1650  PRINT  "CC(RMS)="C1,"PHI="P2*57.296 
1660  Ml=SQR(M3t2+M4t2) 
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1670  P3=ATN(M3/M4) 

1671  X=M3 

1672  Y=M4 

1673  F=P3 

1674  GOSUB  1900 

1675  P3=F 

1680  PRINT  "CM(RMS)="M1,"PHI="P3*57.296 
1690  A1=N1*C0S(A9)*SIN(P1)-C1*SIN(A9)*SIN(P2) 

1700  B1=N1*C0S(A9)*C0S(P1)-C1*SIN(A9)*C0S(P2) 

1710  L9=SQRCAl+2+Bl+2) 

1720  L8=ATN(A1/B1) 

1721  X=A1 

1722  Y=B1 

1723  F=L8 

1724  GOSUB  1900 

1725  L8=F 

1730  A2=N1*SIN(A9)*SIN(P1)+C1*C0S(A9)+SINCP2) 

1740  B2=N1*C0S(P1)*SIN(A9)+C1*C0S(A9)*C0S(P2) 

1750  D9=SQR(A2+2+B2f 2) 

1760  D8=ATN(A2/B2) 

1761  X»A2 

1762  Y=B2 

1763  F=D8 

1764  GOSUB  1900 

1765  D8=F 

1770  A3=M1*SINCP3)-0.75*N1*SIN(P1) 

1780  B3=M1*COS(P3)-0.75*N1*COS(P1) 

1790  M9=SQR(A3f2+B3+2) 

1800  M8=ATN(A3/B3) 

1801  X=A3 

1802  Y=B3 

1803  F=M8 

1804  GOSUB  1900 

1805  M8=F 
1810  PRINT 

1820  PRINT  TAB 10" AERO DYNAMIC  COEFFICIENT  SUMMARY" 
1830  PRINT 

1840  PRINT  "CL="N5"  +"L9"SIN(W1*T+"L8*57.296")" 

1850  PRINT 

1860  PRINT  "CD="C5"+"D9"SIN(W1*T+"D8*57.296")" 

1870  PRINT 

1880  PRINT  "CM(C/4)="M5"+"M9"SIN(W1*T+"M8*57.296")" 
1890  GOTO  90 

1900  IF  X<0  AND  Y<0  THEN  1920 

1910  IF  X>0  AND  Y<0  THEN  1920 

1915  GOTO  1930 

1920  F=F+PI 

1930  RETURN 

1940  STOP 


STEADY  DATA  OUTPUT 


RUN  NO.  IS  123 

COEFFICIENT  BUILDUP  IS 
CNU=  0.3480 
CNL=-0.4281 

TOTAL  CN=- 0.0802 
CCF=-0.0339 
CCR=  0.0721 

TOTAL  CC=  0.0382 
CMU=  0.1092 
CML--0.2840 
CMF-  0.0017 
CMR=  0.0038 

TOTAL  CM=-0.1693  AT  T.E. 


CL(MEAN)=-0.0765  CD(MEAN)=  0.0451 
ALPHA  =-4.9991  DECS 


(CM)C/4 (MEAN) =-0.1092 


UNSTEADY  DATA  OUTPUT 


CN(R^^S)*SIN(PHI)=-l.  538  5 
CC(RMS)*SIN(PHI)=-0.0312 
CM(RMS)*SIN(PHI) =-0.6336 

CN(RMS)*COS (PHI) =-1.4968 
CC (RMS) *COS (PHI) =-0.0042 
CM (RMS) *COS (PHI) =-0.8581 


CN(RMS)=  2.1465 
CC(RMS)=  0.0315 
CM(RMS)=  1.0667 


PHI=225.7868 
PHI=262. 2760 
PHI=216.4451 


AERODYNAMIC  COEFFICIENT  SUMMARY 


CL=-0.0765  + 2.1405  SIN(W1*T+  225.8305  ) 


CD=  0.0451  + 0.1629  SIN(W1*T+  39.2111  ) 

CM(C/4)=-0. 1092  + 0.5836  SIN(W1*T+  63.0435  ) 


1 

RUN  ND.  H230I 

RNBLE  DF  HTTHCK=-5  DEB 

C«c=’0 . 00 
U=100  FT/SEC 

i TAP 

PRESSURE 
♦ CDEFFICIENT 

TAP 

PRESSURE 
* CDEFFICIENT 

1 

0.230 

2B 

-0 . 0B0 

1 ^ 

0.537 

29 

-0. 115 

3 

H 

S 

0.S2B 

30 

-0.0EE 

0.BB7 

31 

-0 . 1 2B 

0.B1E 

32 

-0.107 

E 

0.70E 

33 

-0.133 

7 

0.HH3 

3H 

-0.171 

a 

0.  15E 

35 

-0 . 1 BE 

1 a 

0.00E 

3E 

-0.201 

10 

-0. IB7 

37 

-0.2E1 

1 1 

-0.259 

3B 

-0 . 2B2 

! 12 

-0.353 

39 

-0 . 2BH 

13 

-0 . Han 

M0 

-0 . 339 

IH 

-0.EI 1 

HI 

-0 . 3HB 

15 

-0.G55 

H2 

-0.391 

IE 

-0.E0H 

H3 

-0.503 

17 

-0.592 

HH 

-0 . 503 

IB 

-0.591 

H5 

-0.5B9 

i 19 

i 20 

-0.HEB 

HE 

-0.E72 

-0.HBE 

H7 

-0.B7E 

21 

-0 . 390 

HB 

-1.325 

22 

-0.2MH 

H9 

-1.399 

23 

-0.230 

50 

-0.922 

2H 

-0.075 

51 

-1 .H02 

25 

-0. IHH 

52 

-0 . 9H0 

2E 

-0. IGH 

53 

-0. 195 

27 

-0. I IB 

52 


PRESSURE  CDEFFICIENT  PATH  WITH  STEADY  <U^ 


RUN  ND.  42302 

RNELE  DF  RTTflCK=-5  DEE 


C/4“0. 

li-lZZ 


Z.ZBEE 
FT/5EC 


TAP  ♦ 


PRESSURE 

COEFFICIENT 


S3 


TBP  ♦ 


PRESSURE 

COEFFICIENT 


I 

I .0021 

2B 

0.  1 IE 

2 

Z.974 

29 

-1 . 171 

3 

Z.G39 

30 

-1 .Z72 

4 

31 

-0.EZ4 

5 

Z.233 

32 

-0.379 

E 

Z.Z33 

33 

-0.2E3 

7 

-Z. IBB 

34 

0.013 

B 

-0.44B 

3S 

0.0EE 

9 

“0.E9B 

3B 

0.099 

IZ 

-z.aiB 

37 

0.ZE1 

I I 

-0.9EB 

3B 

0.023 

12 

-1 .074 

39 

-0. 10S 

13 

-1 .20S 

40 

-0.100 

14 

-1.274 

4! 

-0.Z9B 

IE 

-1.343 

42 

-0.172 

IE 

-1 .3SS 

43 

-0.233 

17 

-1 .372 

44 

-0.221 

IB 

-1 .317 

4S 

-0.2E9 

19 

- 1 . 23S 

4B 

-0.21S 

20 

-1 . 190 

47 

-0.279 

21 

-3.371 

4B 

-0.3B3 

22 

-S.SS2 

49 

-0 . 1 77 

23 

-2.B14 

EZ 

-0.020 

24 

-0.4EB 

E! 

0.  1 19 

2S 

-1  . 1S4 

S2 

0.4E3 

2E 

-2.71B 

E3 

0.BI  1 

27 

-2.940 

*• 

1 


■i 
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PRESSURE  CDEFFICIENT  PATH  MITH  5TERDY 

RUN  NO.  MH303 
RNBLE  DF  HTTRCK=-E  DEG 


Cm=0.0GH0 
U=100  FT/SEC 


TAP  * 


PRESSURE 

CDEFFICIENT 


TAP  ♦ 


PRESSURE 

CDEFFICIENT 


I 

1 .000 

2B 

0.339 

2 

0.9B7 

29 

0.203 

3 

0.7B2 

30 

-0.B95 

H 

0.5H1 

31 

-0.279 

5 

0.323 

32 

-0.22E 

G 

0.  15E 

33 

-0. 151 

7 

-0.051 

3M 

0.005 

B 

-0.3B7 

35 

0.01M 

S 

-0.HBH 

3B 

0.023 

10 

-0.E95 

37 

-0.023 

I 1 

-0.B15 

3B 

-0.031 

12 

-0.931 

39 

-0.097 

13 

-1 .0B5 

M0 

-0. 1 IB 

IM 

-1 . IE0 

MI 

-0. 1M0 

15 

-1 .23B 

M2 

-0.22B 

IE 

-1 .23E 

M3 

-0.2EB 

17 

-1 . 1E5 

MM 

-0.311 

IB 

-1 .210 

M5 

-0.29E 

IB 

-1  . 125 

MG 

-0.313 

20 

-1.051 

M7 

-0.3B7 

2! 

-2.5B3 

MB 

-0.5BM 

22 

-H. 1 IH 

M9 

-0.GB7 

23 

-1 .BB0 

50 

-0.231 

2M 

-0.2MB 

5! 

0.037 

25 

-0.71B 

52 

0. 1 17 

2E 

-2.0ME 

53 

0 . 709 

27 

-1 . 151 

54 
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PRESSURE  CnemClENT  DHTfl  WITH  STEADY  U 

RUN  ND.  H230H  g.«0.Z4HB 

HNBLE  DF  RTTRCK=-5  DEE  U=I00  FT/SEC 


PRESSURE  PRESSURE 

TAP  ♦ CDEFFICIENT THP  ♦ CDEFFICIENT 


I 

0.B72 

2B 

0.237 

2 

i .00S 

29 

0.  192 

3 

0.903 

30 

-0.23B 

H 

0.EES 

31 

-0.0B7 

5 

0.53S 

32 

-0. I3E 

E 

0.3I2 

33 

-0. lEH 

7 

0.090 

3H 

-0.0IB 

B 

-0. 197 

35 

-0.03E 

9 

-0.3E2 

3E 

-0.05H 

I0 

-0.SI9 

37 

-0.07E 

I I 

-0.EI9 

3B 

-0. 105 

I2 

-0.775 

39 

-0. 1 19 

I3 

-0.92B 

H0 

-0. 1H3 

IH 

-I .010 

HI 

-0. 1E7 

IS 

-1  .092 

H2 

-0.29H 

IE 

-I .0HE 

H3 

-0.3E2 

I7 

-1 .05H 

HH 

-0.3E7 

IB 

-I .01B 

HS 

-0.3E7 

I9 

-0.BE7 

HE 

-0.3E7 

20 

-0.903 

H7 

-0.53H 

HI 

-1 .9H9 

HB 

-0.777 

22 

-2.99H 

H9 

-0.E13 

23 

-1 .03H 

50 

-0.297 

2H 

-0.0HS 

5! 

-0.379 

HS 

-0.229 

52 

0.090  . 

2E 

-0.7B2 

S3 

0 . 579 

27 

-0.00B 
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PRESSURE  CDEFFICIENT  DHTfl  WITH  STEADY  C/ 


RUN  ND . M2305 

RNELE  nr  RTTRCK=-S  DEB 


PRESSURE 


04=0.0272 
U=100  FT/SEC 


PRESSURE 


0 

[■ 


PRESSURE  CDEFFICIENT  DFITn  WITH  STEADY 

RUN  ND.  M230G  C>*0.0I33 

RNELE  DF  flTTflCK«-S  DEE  U=100  FT/SEC 

PRESSURE  PRESSURE 

TRP  * CDEFFICIENT  TRP  ♦ CDEFFICIENT 


1 

0.S22 

2B 

0.0BG 

2 

0.90B 

29 

0.0B5 

3 

0.9E4 

30 

2.275 

4 

0.B2B 

31 

0.0G3 

5 

0.743 

32 

0.039 

E 

0.S90 

33 

-2.039 

7 

0.319 

34 

-0.0B0 

B 

0.04E 

3S 

-2.150 

9 

-0. 13E 

3B 

-0.220 

10 

-0.31 1 

37 

-0.279 

1 1 

-0.407 

3B 

-0.27B 

12 

-0.4BB 

39 

-0.251 

13 

-0.E04 

42 

-2.291 

IH 

-0.E97 

41 

-2.331 

IS 

-0.7B9 

42 

-0 . 427 

IE 

-0.7G2 

43 

-2.515 

17 

-0.B50 

44 

-0.579 

IB 

-0.Egg 

45 

-2.5E3 

19 

-0.E32 

4E 

-0.E74 

20 

-0.E41 

47 

-2.772 

21 

-0.97S 

4B 

-1 .003 

22 

-1 .309 

49 

-1 . 123 

23 

-0.30E 

50 

-0.454 

24 

0.00E 

51 

-1.275 

2S 

0.0E7 

52 

-2.B3B 

2E 

2.097 

53 

"0 . 290 

27 

Z.0E7 

s: 


RPPENDIX  C 

UN5TERDY  PRESSURE  DRTH 


RNBLE  nr  RTTflCK=“E  DEE  U=100  FT/5EC 


TAP  ♦ 

MERN  PRESSURE 
CDEFFICIENT 

UNSTEADY  PRESSURE 
CDEFFICIENT 

PHRSE 

ANGLE 

1 

0.270 

0.H0H 

1E0 

■ 

2 

0.72E 

0.  IBE 

1E0 

1 

3 

0.90B 

0.073 

230 

1 

H 

0.B12 

0.  lEB 

330 

1 

E 

0.7E0 

0.203 

3HS 

G 

0.B2H 

0.210 

3E3 

7 

0.3B3 

0.22E 

0 

B 

0.  IBl 

0.3BH 

17 

9 

-0.0HB 

0.H00 

23 

10 

-0. IBS 

0.SHE 

32 

i 1 

-0.273 

0.E13 

3H 

12 

-0.3E2 

0.EB1 

37 

13 

-0.HB0 

0.730 

HE 

IH 

-0.EE3 

0.737 

H3 

IE 

-0.E2B 

0.7H3 

El 

IG 

-0.E21 

0.7B0 

E3 

>■. 

17 

-0.E9E 

0.EI3 

EE 

IB 

-0.E97 

0.EG1 

E0 

13 

-0.H9B 

0.E00 

EB 

20 

-0.M21 

0.EB2 

7E 

21 

-0.3E7 

0.HE1 

BB 

22 

-0.313 

0.M00 

103 

23 

-0.2SE 

0.33B 

1 1 1 

2H 

-0.2B3 

0.2EE 

109 

2E 

-Z.27E 

0.  IBE 

32 

2E 

-0.33B 

0.  IBB 

B3 

27 

-0.H2E 

0.  IBB 

BE 

: 

58 

» 

L 

j 

r ^ 


UNSTEADY  PRESSURE  PATH 
C CONTINUED} 

ANBLE  OF  ATTACK=-S  DEB  U-100  FT/SEC 

MEAN  PRESSURE  UNSTEADY  PRESSURE  PHASE 

TAP  ♦ COEFFICIENT COEFFICIENT ANBLE 


2B 

-0.H10 

0.  125 

3H0 

23 

-0.2B5 

0.09H 

295 

30 

-0.23B 

0.  IBl 

290 

31 

-0.257 

0.223 

290 

32 

-0.1B1 

0.31H 

290 

33 

-0. 173 

0.37B 

203 

3M 

-0.I7B 

0.H21 

270 

35 

-0. 192 

0.H7B 

255 

3B 

-0.205 

0.539 

2H0 

37 

-0.22B 

0.B00 

237 

3B 

-0.2H3 

0.5BB 

230 

39 

-0.2B0 

0.BBB 

231 

H0 

-0.2BH 

0.BBH 

227 

HI 

-0.30B 

0.BB0 

22H 

H2 

-0.3BB 

0.B13 

217 

H3 

-0.50B 

0.BI3 

215 

HH 

-0.55B 

0.E00 

210 

H5 

-0.593 

0.530 

205 

HB 

-0.BHB 

0.513 

195 

H7 

-0.B53 

0.HB3 

190 

HB 

-1 .250 

0.H25 

103 

H9 

-1 .300 

0.H05 

175 

50 

-0.9B3 

0.H05 

171 

5! 

-I .H53 

0.H50 

152 

52 

-!  . 125 

0.500 

157 

53 

-0.335 

0.H35 

1H2 
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